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Although nonperturbative functional methods are often associated with low energy 
Quantum Chromodynamics, contemporary studies indicate that they provide reliable 
tools to characterize a much wider spectrum of strongly interacting many-body sys- 
tems. In this review, we aim to provide a modest overview on a few notable applications 
of Dyson-Schwinger equations to QCD and condensed matter physics. After a short 
introduction, we lay out some formal considerations and proceed by addressing the con- 
finement problem. We discuss in some detail the heavy quark limit of Coulomb gauge 
QCD, in particular the simple connection between the nonperturbative Green's func- 
tions of Yang— Mills theory and the confinement potential. Landau gauge results on the 
infrared Yang-Mills propagators are also briefly reviewed. We then focus on less common 
applications, in graphene and high-temperature superconductivity. We discuss recent 
developments, and present theoretical predictions that are supported by experimental 
findings. 
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1. Introduction 

Quite generally, one of the challenges of modern theoretical physics is to explain the 
dynamics of quantum systems. Among them, fermion-fermion interactions are very 
specific and generate a large number of effects — from quark confinement in hadron 
physics, to quark gluon plasma in heavy ion collisions, from superconductivity in 
metals to band structure in crystals. Many different methods have been devised to 
tackle the many-body problem, under a variety of circumstances and at different 
energy scales. In high energy QCD, where asymptotic freedom guarantees that the 
coupling between quarks and gluons can be treated as a small parameter, pertur- 
bation theory has been successfully applied. At intermediate and low momenta, 
different methods are required since the coupling increases as we go to lower en- 
ergies and perturbation theory alone cannot give a good description of the theory. 



1 



25, 2013 2:8 WSPC/INSTRUCTION FILE mpla-rev 



2 Carina Popovici 

One possibility is to employ continuum functi ona l techniques such as renormal- 
ization group^ and Dyson-Schwinger equations 12121 Alternatively, numerical lattice 
simulations represent a potentially exact method, howev er e xtrapolation to chiral 
and infinite volume limits is often a numerical challengePEI The situation is very 
similar in condens ed mat ter electronic systems, where an equivalent program has 
been implemented ISiHSlSI 

This review is organized as follows. In section [2] we briefly review the general 
features of Dyson-Schwinger equations, and proceed in section |3] by discussing the 
confinement problem, with particular emphasis on the heavy quark limit of Coulomb 
gauge QCD. A brief overview of the progress made on infrared propagators of Lan- 
dau gauge Yang-Mills theory is also attempted. Phenomenological findings are not 
discussed here, since they have been presented in many other works, see for example 
Ref. llOl for a recent review. In section|4]we turn to planar condensed matter systems. 
We exemplify the Dyson-Schwinger approach to strongly interacting many electron 
systems with two case studies, graphene and high temperature superconductors 
(HTSs). Our choice is motivated by the fact that these systems share a common 
aspect, namely, they can be described by effective quantum field theoretical models 
that inherit features of both (2+1) dimensional Quantum Electrodynamics (QED3) 
and its four dimensional counterpart. We present results obtained from the Dyson- 
Schwinger equations, and where available the supporting experimental evidences. At 
the more formal level, it is also worth mentioning the similarities between the prop- 
agators of Yang-Mills theory and QED3, i.e., the possibility of power law behavior 
in the deep infrared. A short summary is presented in section [S] 

2. Dyson Schwinger formalism 

Dyson-Schwinger equations built up an infinite set of coupled non-linear integral 
equations that relate the various Green's functions of a quantum field theory. From 
a practical point of view, defining a tractable problem requires the ability to reduce 
these equations to a closed subset, e.g. by making controlled ansatze for the higher 
order Green's functions. In general the resulting coupled integral equations are 
solved numerically, nevertheless it is possible that in particular setups, such as 
heavy quark sector of Coulomb gauge QCD (under certain truncations), or infrared 
limit of Landau gauge Yang-Mills theory, exact analytic solutions become available. 
These aspects will be discussed in more detail in the next section. 

The derivation of the Dyson-Schwinger equations rests on the idea that the 
generating functional of a theory 




(2.1) 



is invariant under the variation of a generic field $„: 




(2.2) 
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Fig. 1. Dyson-Schwinger equation for a generic fermion propagator in a gauge theory. Filled 
circles denote dressed propagators, the empty circle denotes the dressed fermion- gauge boson 
vertex, the spring denotes the gauge propagator and solid lines stand for the fermion propagator. 

where S is the action of the problem at hand, Ss the corresponding source term and 
2?$ denotes the integration over ah fields. To ensure the validity of this equation, 
translational invariance of the measure 2?$ is assumed. From Eq. (|2.2|) . one can 
derive the Dyson-Schwinger equation for any n-point Green's function, by taking 
further functional derivatives with respect to the fields and omitting those terms 
which vanish when the sources are set to zero. 

Among the most investigated Dyson-Schwinger equations is the fermion gap 
equation, which, regardless of the underlying gauge theory, has the diagrammatical 
representation shown in Fig. [T] It essentially states that the inverse of the nonper- 
turbative fermion propagator equals the sum between the inverse of the bare fermion 
propagator and the fermion self energy, which in turn contains dressed fermion and 
gauge boson propagators, and one bare and one dressed fermion-gauge boson ver- 
tex. Formally, the fermion gap equation appears in systems as discrepant as QCD, 
graphene and superconductors, whereby the raw form Fig. [T]has to be customized 
to the particular theory under consideration. 

In addition to formal analogies, there are also important qualitative similarities 
between the systems that we are going to examine. An example is the possibility 
of scaling laws for the propagators in the deep infrared — a feature shared by both 
Landau gauge Yang-Mills theory and QED3 (including high temperature supercon- 
ductors, which can be described with a model that has properties reminiscent of 

QED3).i2inini 

3. Confining solutions from the Dyson Schwinger equations 

Confinement, i.e., the experimental evidence that quarks and gluons must be con- 
fined into color singlet hadronic states, represents, along with chiral symmetry 
breaking, one of the theoretical pillars of low energy hadron physics. A rigorous 
understanding of this phenomenon, i.e., identifying the specific mechanisms that 
act at the level of the underlying (gauge dependent) Green's functions, is to date 
still an open problem. 

As is well known, in order to overcome the difficulties related to the definition 
of the generating functional Eq. (|2.1|) . the QCD action needs to be supplemented 
by a gauge fixing term, along with the associated ghost fields. Among various op- 



25, 2013 2:8 WSPC/INSTRUCTION FILE mpla-rev 



4 Carina Popovici 

tions, Coulomb gauge diAf — has properties that quahfy it as an efficient choice 
to study nonperturbative phenomena: in this gauge, Gauss law is naturally built 
in such that gauge i nvarian ce is fully accounted for, the total color charge is con- 
served and vanishing j 1 3| 1 4| most importantly, a natural picture of confinement 
emerges. Furthermore, Coulomb gauge allows for different manifestations of con- 
finement depending on the specific formalism, as shall be discussed below. 

In the Hamiltonian formalism. Coulomb gauge is implemented simultaneously 
with Weyl gauge Aq — ^.^'^ ^'^ In this approach, the Dyson-Schwinger 

equations for the ghost and spatial gluon propagators have been s olved, analytica lly 
in the infrared and numerically in the whole momentum regime I f6 | 17 | 20 | 21 | 22 | 
Ref. [23] recent re sults at fi nite temperature are presented). Similar to Landau gauge 
scaling solutions ,1^312326] ^j^g gj-^ost and gluon dressing functions exhibit power law 
behaviors, {p^Y^'^ and [p^Yn' . A general infrared power law analysis has yielded 
the following relation between the infrared exponentj^ (with Landau and Coulomb 
gauge differing only in the number of dimensions d) ' ^'^ I ^O i 

4 - d 

5gi + 25gh = ^ ■ (3.3) 

In the case d = 3, which corresponds to equal-time Green's functions in Coulomb 
gau ge, two power law exponents k := —5gh have been found, k*^^^ = 0.398 and kP'^ = 
1/2 EQl. The energetically favored solution is the most singular, which produces a 
ghost propagator dressing function diverging as l/|fc|. In turn, this gives rise to a 
strictly linearly rising static quark potential at large distances.! ^ The case d = 4 
trivially leads to the Landau gauge scaling relation 5gh — —6gi/2^^^ (see also the 
discussion at the end of this section). 

In the functional formalism (as applied to Coulomb gauge QCD), apart from 
the ghost and transversal spatial gluon propagator, there is a third propagator, 
the temporal gluon propagator, which plays an important role in uncovering the 
confining potential between (static) quarks. To see how this comes about, let us 
first briefly review the main characteristics of the functional formalism. 

In general, functional Dyson-Schwinger studies in Coulomb gauge are plagued by 
the so-called energy divergence problem, i.e., the unregulated divergences generated 
by ghost loops, as in the one-loop integral / dk^ J d^k[{k — The usual 

dimensional regularization fails in this case, although the full set of such loops 
should cancel in the Dyson-Schwinger equations. While such cancellations have 
been isolated up to two loops in perturbation theory,l2Sl29l30l3T] ^j-^gy exceedingly 
difficult to pin down in the full tower of Dyson-Schwinger equations. This problem 
can be bypassed by converting to first order functional formalism, i.e., linearizing 
the chromoelectric term in the action via the auxiliary field 7? 



exp <{ i / d'^x^E-E \ = / Vncxp \i j d'^x 



TT-TT — TT-E 

2 



(3.4) 



''Note, however, that imposing Weyl gauge leads to the fact that Gauss law is no longer arising 
from the equations of motion, but has to be imposed as an external constraint. 
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Fig. 2. Homogeneous Bethe-Salpeter equation for quark-antiquark bound states. Solid lines rep- 
resent the (dressed) quark propagators, dashed hnes denote the corresponding bound state, and 
blobs denote the Bethe— Salpeter vertex function. The box represents the Bethe-Salpeter kernel 
which in the heavy quark limit reduces to the exchange of a single temporal gluon. 



The TT field is subsequently split into transverse and longitudinal components (details 
are provided in Ref. [5^ . With this new fields, the QCD action is rewritten such 
that the ghost field, i.e., the Faddeev- Popov determinant, cancels against the inverse 
functional determinant that stems from resolving the chromodynamical equivalent 
of the Gauss law. The resulting action contains only the 'would-be-physical' degrees 
of freedom, the transverse A and 7? fields (which classically would correspond to the 
configuration variables and their momentum conjugates), whereas the unphysical 
ghosts are formally eliminated. The shortcoming, however, is th at we end up with a 
nonlocal Lagrangian that is unsuited for practical applications .tl^ISll Nevertheless, 
the nonlocal formulation can provide important guidelines to the local theory. 

It is important to realize that identifying the mechanisms that lead to the can- 
cellations of the unphysical components is essential for a thorough understanding of 
the theory. For example, the (unphysical) energy-independent ghost loop mentioned 
above should be canceled by the temporal component of the gluon propagator, and 
this implies that the temporal gluon itself must have an energy-independent part.^^ 
This argument is supported by lattice results and analytic findings In the 
limit of the heavy quark mass, this information has been used as input to derive 
a relation between the temporal gluon propaga tor and the nonperturbative scale 
associated with confinement (the string tension). We will return to this relation 
below. 

Explicitly, heavy quark limit refers in this case to a heavy quark mass expansion 
of the QCD action at leading orderJB with the truncation of the Yang-Mills sector 
to include only dressed propagatorsJB In this framework, a full nonperturbative 
study of the gap and Bethe-Salpeter equations has been performed. Solving the 
gap equatioio (generically represented in Fig. [T|) yields the following solution for 



The QCD action is expanded in powers of the inverse quark mass by means of a heavy quark 
transformation adapted from Ref. 1371 

"^A similar Coulomb gauge truncation scheme with arbitrary quark mass h as also been studied, 
and which explicitly reproduces the heavy quark limit presented here.EHESl 

"^The full gap equation within Coulomb gauge first order formalism (without the mass expansion) 
has been derived in Ref. 1401 
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the heavy quark propagator: 

W^giko) = jr ~V ^ - r (3.5) 

[ka - m-Ir + ie\ 

where m denotes the quark mass and Ij. is an (imphcitly regularized) constant. 
Since this propagator has a single pole in the complex fco plane, it follows that the 
closed quark loops are suppressed, and this implies that the theory is quenched in 
the heavy mass limit: 

dk[)— — — — ~ 0. (3.6) 

[fco - m - Ir + lej [ko+Po - "m- X,. + ie\ 

By using the associated Slavnov-Taylor identity, it is easy to show that the temporal 
quark-gluon vertex remains nonperturbatively bare, TqqAo — ^fqAo — 5^°: with 
being the hermitian generators of the gauge group. Moreover, the kernel of the 
Bethe-Salpeter equation, depicted in Fig. [21 reduces to a single gluon exchange: 

K{k)^^,AoWS^{k)T^,lAo^ (3.7) 

where Woo(fc) is the temporal gluon propagator (a complete description and deriva- 
tion is provided in Ref. [36|). With the kernel Eq. (|3.7p . the bound state energy 
between a quark and an antiquark at leading order in the mass expansion reads: 

^g' (^^o°('^) - ^A/e'"-^l , (3.8) 

where the index V denotes an implicitly regularized integral, Cf is the Casimir 
factor. Cm denotes an (unknown) color factor assigned to the quark-meson vertex, 
and \x\ represents the separation between the quark and the antiquark. To identify 
Cm, we recall that since a quark cannot live as asymptotic state,fI3l the bound 
state energy Pqq must be linearly rising, or otherwise the energy is infinite and the 
quark-antiquark system is not allowed. With a temporal gluon dressing function 
more divergent than l/|w|, it follows that Cm = Cp and hence a finite solution of 
the Bethe-Salpeter equation exists only for color singlet states. Assuming that in 
the infrared 

1^00 (c3) = 4' (3-9) 

wh ere X is a combination of constants (and g^X is a renormalization group invariant 
Eq. (I3H) reduces to 

Pqq:^a\x\^^^^^\x\- (3.10) 
<37r 

This result provides the link (at least at leading order) between the physical string 
tension a and the nonperturbative Yang-Mills sector of QCD. It corresponds to the 
only physical solution (a color singlet bound state of a quark and an antiquark) or 
else the energy of the system is divergent. A similar calculation performed for the 
diquark Bethe-Salpeter equation has shown that diquarks are confined for Nc — 2 
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K(d: 




Fig. 3. Faddeev equation for three-quark bound states. The box depicts in this case the diquark 
kernel, the elhpse denotes the Faddeev vertex function, and the triple lines represent the three- 
quark bound state. Same conventions as in Fig. [2] apply. 



colors, which corresponds to a confined, antisymmetric bound state of two quarfcs 
(the SU{2) baryon), and otherwise there are no physical states. Finally, a thorough 
investigation of the Dyson-Schwinger equation for the full nonperturbative four- 
point quark-antiquark Green's functions has evidenced the separation of the physical 
and unphysical singularities of the Green's function, and that the physic al p ole 
coincides with the pole Eq. (|3.8[) of the homogenous Bethe-Salpeter equation.!^ An 
analogous investigation for the diquark system has shown that the resonant pole of 
the four-point Green's function is attached to the only physical solution, for Nc = 2 
colors, corresponding to color antisymmetric and flavor symmetric configuration. 

Following the same approach, a similar analysis has been carried out for three 
quark systems by means of Faddeev equation, depicted in Fig. [3p21 Since in the 
quark-diquark model the binding energy between quarks is assumed to be provided 
essentially by two-quark correlations, in the Faddeev equation only the permuted 
two-quark kernels are employed, whereas the three-body irreducible diagrams 
are neglected. Just like for qq states, the diquark kernel reduces to the ladder ap- 
proximation. 



K,. 



^(rf)(fc) = r|,AoW^oo(^)q,Ao- (3-11) 

By working in the symmetric case (equal spatial separations between quarks), the 
following three-quark bound state energy has been determined as a solution of the 
Faddeev equation 



P 



3q 



3m + - g 
2^ 



did 

(2^ 



(3.12) 



where Cb denotes the color factor attached to the quark-baryon vertex. As before, 
the only possibilities are either that the 3q system is confined, or it is physically not 
allowed. An infrared confining solution requires that the condition Cb = Cf/2 is 
satisfied, provided that the temporal gluon propagator is infrared enhanced. This 
implies that the baryon is a color singlet bound state of three quarks and otherwise 
the energy of the system is infinite. Inserting Cb — C'f/2, along with the infrared 
temporal gluon propagator Eq. p.9|) . into Eq. p.l2|) . yields the following result for 
the bound state energy of the three quark system: 



3q 



^ig^CpX 



(3.13) 
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As in the case of qq systems, the 'string tension', i.e., the coefficient of the three- 
body hnear confinement term, is directly related to the nonperturbative temporal 
gluon propagator. Comparing the above result with Eq. p.lO|) . it appears that the 
string tension corresponding to the qqq system is 3/2 times that of the qq system. 
The 3m term is a reminiscence of the heavy quark transformation which enters the 
Faddeev equation via the heavy quark propagator, Eq. p.Sp . 

Let us now turn our attention towards Landau gauge 9^^^ = 0. Unlike Coulomb 
gauge, where one is faced with severe technical difficulties introduced by noncovari- 
ance, Landau gauge remains popular due to its covariance, but also because of cer- 
tain technic al si mplifications (the ghost-gluon vertex remains approximately bare in 
the infrared this gauge, the coupled system of Dyson-Schwinger equations 

for ghost and gluon propagators has been solved, see also the discussion above 
Eq. (|3.3p . It was found that the gluon propagator vanishes at zero momentum, 
whereas the ghost propagator is infrared enhanced, in agreement with the ghost 
dominance picture of Gribov and Zwanziger,^*^ '^'^ ^"^ and Kugo and Ojima.^^ '^^^ 
With the truncation to employ a bare ghost-gluon vertex, the scaling coefficient 
K « 0.595 has been calculated independently in Refs. [27l and [26pl Bv combining in- 
formations stemming from both Dyson-Schwinger and renormalization group equa- 
tions, it has been shown that the Landau gauge scaling solution is unique (however, 
this should not be confused with the uniqueness of the solution itself).!^ 

On the lattice, it eventually become clear that the gluon propagator is not van- 
ishing at small momentum; instead, a massive behavior has been seen, although 
some papers were still in agreement with a vanishing gluon propagator, see for 
example Rcf. 52 and references therein. These findings support the so-called decou- 
pling solution — a second solution derived with Dyson-Schwinger methods, where 
the gluon and ghost propagate quite differently: at small momenta, the gluon prop- 
agat or become s finite instead of going to zero, whereas the ghost dressing remains 
finite.lSniMISSI In order to understand the connection between the two solutions, it 
is necessary to examine the renormalization of the ghost propagator. A finite ghost 
dressing function at zero momentum generates a continuous set of decoupling so- 
lutions, whereas if an infini te g host dressing at zero momentum is employed, the 
scaling solution is recovered.'^' Finally, we note that a numerical analysis of gluon 
an d gh ost propagators in the complex momentum plane has been recently carried 



out.E3 

Importantly, both the decoupling and scaling solutions fulfill the criterion for a 
confining Polyakov loop potential, meaning that quarks are confined.^^ This crite- 
rion depends only on the asymptotic part, and as it turns out, in the actual calcula- 
tions the dynamics of the system is driven by the non-perturbative mid-momentum 
regime (where both solutions agree), and not the deep infrared. 



^In fact, it turns out that this result changes very little if the vertex is dressed in the infrared.Li 
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4. Manifestations of (2+1) dimensional QED in condensed matter 
physics 

Although when studying gauge theories it seems natural to consider a theory set 
in 3 + 1 dimensions, it is nevertheless possible that some insight might be gained 
by restricting to 2 + 1 dimensions. An example is QED3, which apart from having 
intrinsically interesting features, such as superrenormalizability,^'^ is also attractive 
for its ability to describe the behavior of certain materials that lay at the border 
between particle and condensed matter physics, such as graphene and high temper- 
ature superconductors. The quantum field theoretical modeling of these materials 
has provided a description of phenomena that are nonperturbative in essence, e.g. 
dynamical mass generation via chiral symmetry breaking. 

Before turning our attention to condensed matter systems, we note that the 
infrared behavior of Landau gauge QED3 propagators is given by power laws, just 
like the scaling solutions of Landau gauge Yang-Mills theory (see also the discus- 
sion from the previous section). The QED3 power laws have been derived from an 
infrared analysis of the coupled system of fermion and photon Dyson-Schwinger 
equations 



4.1. Graphene 

Graph ene is a monatomic layer of carbon atoms arranged on a honeycomb 
lattice.^ ^l^-*^ ! ^^! Its rem arkable electro nic pro perties, such as unconventional quan- 
tum Hall effect,IS2I^ Klein tunneling I65|66| charge confinement,!^ qualify it as 
one of the most promising materials for future nanoscale devices. Its synthesis in 
2004'^^'^^' was awarded the Nobel Prize in physics in 2010. 

The fundamentally new characteristics of this system are theoretically explicable 
in terms of the low energy behavior of the charge carriers: instead of the familiar 
Schrodinger equation, the electrons in gra phene are described by the relativistic 
Dirac equation for massless fermions.^^^I^ The resulting single particle dispersion 
relation is linear in the momentum, Ek = ±VF\k\, where vp is the Fermi velocity 
(about 300 times smaller than the speed of light) and k — [k^, ky) is the fermion 
momentum, measured relative to the inequivalent corners of the Brillouin zone. 

From a technological point of view, graphene based nanoscale devices require the 
ability to control the so-called charge confinement, i.e., the clustering of electronic 
charge such that the Klein effect can be overcome.^^ Theoretically, a low energy 
effective model that is well suited to describe variou s disorder phenomena has been 
constructed by means of a chiral gauge theory.f^2l23| ^jjjg ]^a,s been supplemented by 

scalar and gauge fields that account for the self-interaction of the carbon background 

74. 

and the mean self-interaction of the Dirac fermions:-*-* In this framework, charge 
confinement has been described by modeling one-dimensional defects, associated 
e.g. with chemical bonding of foreign atoms to carbon atoms, as potential barriers 
which break sublattice symmetry.^^ 

The single-particle picture, however, changes quite dramatically if the long range 
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Coulomb interaction is taken into account. Since fermions are confined to live in the 
plane, whereas the gauge bosons live in three spatial dimensions, one is confronted 
with the problem of coupling the two-dimensional current with a three-dimensional 
gauge potential. This difhculty is solved by integrating the usual photon propagator 
over the third momentum dimension (known as dimensional reduction), which leads 
to an atypical behavior of the gauge propagator, namely, it goes like 1/fc, instead of 
This anomalous gauge propagator induces the renormalization of the electron 
self-energy, which is not present in ordinary QED3. 

Various theoretical studies have predicted that Coulomb interaction may gener- 
ate dynamically a finite mass gap, such that graphene und ergoes a phase transition 
from semimetal to insulator above some critical coupling q;^I 76 | 77 | 78.8 , 79 ^j-^is possi- 
bility is motivated by the large value of the 'bare' coupling a = e^/ehvp of graphene 
on a substrate with dielectric constant e (for suspended graphene ao = 2.19). 

On the other hand, most recent experiments find no signature of an insulating 

on 01 

phase in suspended graphene. Measurements indicate, however, that the real 

dispersion relation is logarithmic, instead of linear, near the neutrality point. This 

reshaping of the Dirac cones is caused by the charge-carrier density dependent 

renormalization of the Fermi velocity (induced by the Coulomb interaction). In fact, 

the running of the fermion velocity with the energy was already predicted in earlier 

82 

renormalization group calculations," where it was found that the Fermi velocity 
grows logarithmically without bound, until retardation effects become important 
enough to invalidate the instantaneous Coulomb approximation. 

To analyze the problem of the gap generation, it is useful to inspect the Dyson- 
Schwinger equation for the fermion propagator Fig. [T] Following the usual proce- 
dure, this equation is converted into a system of coupled integral equations for the 
Fermi velocity dressing function A(p) and the gap function A(p), which characterize 
the fermion propagator. 

Sip)-' = 7% - vfA{p) -f'p, - A{p). (4.14) 

At the critical coupling etc, where the phase transition should take place, the 
system is well described by bifurcation theory. Applying bifurcation theory amounts 
to linearizing the equations around the critical point, i.e., neglect quadratic and 
higher order terms. By using a one-loop photon propagator, it has been found that 
the function A{p) is logarithmically divergent ,lS3]jjj agreement with the experimental 
evidences Explicitly, it reads (g = aNfir/A) 

A{p) = 1 + {- - 2g -f c{g){g' ~ 1)] In | + f{g)] , (4.15) 

where c{g) and f{g) are functions of the coupling, Nf ~ 2 for monolayer graphene, 
and the physical cutoff A is determined by the size of the graphene's Brillouin zone. 
This result is then used in the bifurcation analysis of the gap equation, which around 
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the critical point reduces to: 

[ <Pk 1 A(fc) J k ,,,, \ 
vpe J {2TTy k\p-k\A[k) \\p-k\ J 

In the above, J is a piecewise function that includes the effects of the Fermi velocity 
dressing function. Importantly, for determining the critical coupling corresponding 
to the semimetal-insulator transition, not only the logarithmic but also the regu- 
lar contributions to the function A{p), i.e. the function f{g), become significant.'^ 
For the critical coupling it is obtained ac — 2.85, which is larger than the bare 
coupling of suspended graphene ao = 2.19. These findings are in agreement with 
the experimental observations that suspended graphene does not undergo a phase 
transition. Hence, it appears that the logarithmic renormalization of Fermi ve- 
locity strongly weakens the Coulomb interaction, favoring the persistence of the 
semimetal phase in suspended graphene.!^ 

Despite the encouraging results, a complete understanding of the physical pic- 
ture requires a full nonperturbative investigation, by including the Dyson-Schwinger 
equation for the gauge propagator. Since nonperturbative contributions to both pho- 
ton propagator and vertex functions have led to significant corrections in ordinary 
QED3, it is likely that including these effects will lead to further corrections of the 
critical coupling. Furthermore, a proper comparison of our value for the critical 
coupling with lattice calculations is impeded by the details of the cutoff procedure 
in the ultraviolet, as lattice simulations have been performed on both 'standard' 
squared^ and hexagonal (physical) lattices. Finally, as well known from ordi- 
nary QED3, finite volume effects may also play a role, see for example Ref. 85 and 
references therein. 

When comparing to the experiment, one has to keep in mind that other type of 
terms, such as four-fermion contact interactions, may also become important and 
thus have to be included in the model Lagrangian, in addition to the long-range 
Coulomb interaction. First studies in this direction have already been undertaken 
in Refs. 76^ ^86* 3!L Another important line of study is to explore the effects of a 
finite chemical potential on the Fermi velocity — experimentally this corresponds 
to a finite density o f ch arge carriers, i.e., chemical doping. This problem is currently 
under investigation.!^ 



4.2. High-temperature superconductivity 

Due to their properties apparently incompatible with conventional metal physics, 
high-Tc cuprate superconductors can be regarded as another example of QED3 'in 
the lab'. A particularly striking feature is the pseudogap — a third phase that lays 
at the boundary between the nonsuperconducting (in this case, ass ocia ted with the 
antiferromagnetic spin density wave) and superconducting phases.!^ The physics 
of the pseudogap phase is dominated by strong pairing fluctuations, as has been 
experimentally demonstrated by measuring the so-called Nernst signal.!^ 
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Describing the pseudogapped phase theoreticahy is a chaUenge due to the com- 
pheations introduced by the scattering of the fermions off the vortices (topological 
defects), which are otherwise bounded in the superconducting phase!^ A way to 
circumvent this problem has been put forward by Franz and Tesanovic, who adopted 
an 'inverted' paradigm and proposed to start from the 'conventional' superconduct- 
ing p has e and try to find the mechanism that describes the transition to the normal 
stateP2l In this approach, the fermion field is gauge transformed to a new field 
which the authors call 'topological fermions'. The pseudogap phase of the topo- 
logical fermions is then associated with the chiral symmetric QED3, whereas the 
transition to the antiferromagnetic phase is indicated by the generation of a dy- 
namical mass for the guasiparticles.™™ The generation of such a mass term 
results from interactions of the fermionic quasiparticles with topological excitations 
described by the gauge fields of QED3. 

An open question is whether at zero temperature a pseudogap associated to the 
onset of antiferromagnetic spin density wave instability is formed in the first place, 
or if the system is driven directly into the antiferromagnetic phase. The answer 
lays in the critical number of 'flavors' of fermions where the system undergoes 
a phase transition from the chirally broken into the chirally symmetric phase. This 
has to be compared with the physical number of flavors Nf = 2, which denotes 
the number of low energy quasiparticles located at the four nodal points of the gap 
function, see for example Ref. for details. If < 2, the theory will flrst go 
through the pseudogap phase; if, on the other hand, Nj > 2, the system is driven 
directly from the superconducting into the antiferromagnetic phase upon further 
under doping. 

Furthermore, since experimental observations suggest that most of the HTS ma- 
terials are strongly anisotropic due to the large difference between the Fermi velocity 
and a seco nd velo city va (related to the amplitude of the superconducting order 
parameter) j ^^ l ^^ l a complete theoretical description requires taking into account this 
intrinsic anisotropy as well. At the level of the Lagrangian, including the anisotropy 
amounts to implementing a nontrivial metric that collects the anisotropic fermion 
velocities 

Much work has been devoted to determining the critical number of flavors at 
which chiral symmetry is broken, in the isotropic l^7 | i^8|99 111J ^ ^^^1 as anisotropic 
| 93 | 89 | 100 112] versions of QED3. Studies carried out for isotropic QED3 indicate that 
N'^ lays around 3.5 — 4. In the following, we concentrate on anisotropic QED3, 
and briefly review the results obtained within the Dyson-Schwinger formalism. In 
a first study, the Dyson-Schwinger equation for the fermion propagator has been 

no 

solved in the limit of small anisotropic velocities. It has been found that in this 
case there should be essentially no difference between the universal behavior of 
isotropic and anisotropic QED3. Later on, a more sophisticated fermion dressing 
has been i ncluded , and again only the case of small anisotropic velocities has been 
considered.l^SlinQ] jjj fj^^ large- A^/ limit, it has been found that the critical number 
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of flavors where th e the ory suffers the chiral phase transition stays the same. In 
a subsequent studyji^ a quantitative criterion that relates the critical number of 
fermion flavors and the strength of the photon-fermion coupling has been proposed, 
namely, it was conjectured that the dynamical fermion mass will be generated as 
soon the gauge field strength is larger than some threshold value which can be 
determined from the S-matrix for fermion-fermion scattering. It was concluded that 
velocity anisotropy does affect the number of critical fermion flavors at which chiral 
symmetry is broken due to mass generation. A recent reexamination of the situation 
going beyond the small anisotropy expansion, and using a power law ansatz for the 
photon propagator, ha s shown sizable deviations of the critic al numbe r of flavors 
from the isotropic case,^i^ll^ in agreement with lattice results. H^^HDH 

5. Summary 

In this brief review we have presented a few selected applications of Dyson- 
Schwinger approach to strongly interacting fermion systems. In the first part we 
have considered the confinement problem, concentrating on the heavy quark limit of 
Coulomb gauge QCD. With the truncation to include only dressed Yang-Mills prop- 
agators, we have sketched the derivation of a simple relation between the Green's 
functions of Coulomb gauge Yang-Mills theory and the quark confinement poten- 
tial, for quark- antiquark and three-quark systems. Confining (finite energy) solu- 
tions exist only for color singlet meson/baryon states, and otherwise the systems 
have infinite energy. Further, we have seen that in the heavy quark limit, the rain- 
bow approximation to the quark gap equation and the ladder approximation to the 
associated Bethe-Salpeter and Faddeev equations remain nonperturbatively exact. 
In the remainder of this part we have briefiy reviewed Landau gauge results on the 
infrared Yang-Mills propagators. 

In the second part we have moved on to less commonplace applications, in 
graphene and high temperature superconductors. These planar condensed matter 
systems are well described by effective QEDs-like theories, which are tailored to 
accommodate the particularities of the system under consideration. We have re- 
viewed selected topics studied with Dyson-Schwinger methods, which include the 
problem of dynamical gap generation by long-range Coulomb interactions in sus- 
pended graphene, and the pseudogap phase in high temperature superconductors. 
In graphene, including renormalization effects on the Fermi velocity has lead to the 
conclusion that the Coulomb interaction is strongly weakened near the charge neu- 
trality point, preventing the emergence of a gapped phase. These results agree with 
the experimental findings which indicate the persistence of the semimetal phase. In 
high-temperature superconductivity, where the physics of the pseudogap is strongly 
dependent on the number of fermion species, it has been shown that the inherent 
spatial anisotropy leads to a significant alteration of the critical number of fla- 
vors where the phase transition from chirally symmetric to chirally broken phase 
takes place (as compared to isotropic QED3). Least but not last, it is important to 



25, 2013 2:8 WSPC/INSTRUCTION FILE mpla-rev 



14 Carina Popovici 



appreciate the possible similarities between different theories, i.e., the scahng solu- 
tions ol Landau gauge propagators of Yang-Mills theory and QED3 . In the context 
of high-Tc superconductivity, a power law ansatz for the photon propagator has 
been successfully employed to study the dynamical breaking of chiral symmetry. 
Whether the gauge field in graphene is characterized by a similar behavior remains 
to be clarified in the future. 
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